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Multivariate fractional Poisson processes and compound sums 


Luisa Beghin* Claudio MaccP 


Abstract 

In this paper we present multivariate space-time fractional Poisson processes by considering 
common random time-changes of a (finite-dimensional) vector of independent classical (non- 
fractional) Poisson processes. In some cases we also consider compound processes. We obtain 
some equations in terms of some suitable fractional derivatives and fractional difference opera¬ 
tors, which provides the extension of known equations for the univariate processes. 
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1 Introduction 


Typically fractional processes are defined by considering some known equations in terms of suitable 
fractional derivatives. In this paper we deal with fractional Poisson processes which are the main 
examples among counting processes; here we recall the references hd, m, 0 , 0 , m and m 
(we also cite [TO] and m where their representation in terms of randomly time-changed and 
subordinated processes was studied in detail). Moreover, as pointed out in [20], a class of these 
processes demonstrate the phenomenon of anomalous diffusion (i.e. the variances of the process 
increase in time according to a power f 7 , with 7 yf 1); this aspect was also highlighted in [6] where 
the authors refer to the long-range dependence property (they also present some applications in ruin 
theory where the surplus process of an insurance company is modeled by a compound fractional 
Poisson process). 

The aim of this paper is to present m -variate space-time fractional (possibly compound) Poisson 
processes; in this way we generalize some results in the literature for univariate processes, which 
can be recovered by setting rn = 1. Often closed formulas for fractional Poisson processes are given 
in terms of the Mittag-Leffler function, i.e. 


E a A x ) : = X 

r>0 


T(ar + P) 


(1) 


(see e.g. [18], page 17). 

We start with the simplest case, i.e. the multivariate version of the space-time fractional 
Poisson process in [TO]. In particular we consider the time-change approach in terms of the stable 
subordinator and of its inverse (see (3.18), together with (3.1), in [2]; see also [22]). So we introduce 
the following notation: for u £ (0,1), let { A v (t ) : t > 0} be the stable subordinator and let 
: t > 0} be its inverse, i.e. 


C u (t) := infjz > 0 : A u (z) > t}. 
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In what follows we denote the continuous density of L u {t) by and the continuous density of 

A u (t) by Stable subordinators are well studied in the references on Levy processes (see e.g. 

P 0 and [H]) ; for the inverse of stable subordinators, we recall 0 , m and 13- 

Definition 1.1 Let {{ Ni(t ) : t > 0} : i £ {1,... , m}} be rn independent Poisson processes with 
intensities X\,... , X m > 0, respectively, and set 

N (£) := (t),..., N m (t)). 

Then, for r),v £ (0,1]. we consider the m-variate process {N ri ’ u {t) : t > 0} defined by 

N^ v (t) := N (A v (L u (t))), 

where { N(t ) : t > 0} , {A v (t) : t > 0} and {^(t) : t > 0} are three independent processes. Moreover 
we also consider the cases rj = 1 and/or v = 1 by setting A l {t) = t and L 1 (t) = t, respectively; 
thus, in particular, {IV 1,1 (£) : t > 0} coincides with { N(t ) : t > 0}. 

We remark that {{N/ U (t) : t > 0} : i £ {1,..., m}} in Definition 11.11 are conditionally indepen¬ 
dent given { A ri (L L '(t)) : t > 0} (except for the case rj = v = 1 where they are independent). 

Throughout this paper we deal with m-variate processes and we use the notation a = (ai,..., a m ) 
for ?n-dimensional vectors. For instance we often write k > 0 where k\,... ,k m are nonnega¬ 
tive integers (because we deal with processes with nonnegative integer-valued components) and 
0 = (0,..., 0) is the null vector. Moreover we write: a < b (or a > b) to mean that 
(or ai > bi ) for all i £ {1 ,..., m}; a -< b (or a >- b) to mean that a* < 6, (or a* > bfi for all 
i £ {1,... ,m}, but a / b. Finally we remark that the probability generating functions assume 
finite values when their arguments u belong to [0, l] m but, in some cases, the condition u £ [0, l] m 
can be neglected or weakened (for instance, when rj = 1, this happens for the probability generating 
functions in (plj and ([5]); in the first case the finiteness of G'i(tti),..., G m {u m ) is also needed). 

Our results mainly concern the state probabilities {{p r / u {t) : k > 0} : t > 0} defined by 

Ph’ u (t) := P(N v ’ u (t ) = k ) for all integer k \,..., k m > 0. (2) 

We also consider two generalizations of the process {N T1,u (t) : t > 0} in Definition 11.11 we 
mean the multivariate space-time fractional compound Poisson process (see Definition 11.21) and the 
multivariate version of the process in m, where we have a general subordinator associated to a 
Bernstein function / in place of the stable subordinator {A v (t) : t > 0} (see Definition 11.31) . We 
start with the first one. 

Definition 1.2 For ipv £ (0,1], let {C v,u (t) :t > 0} be defined by 

C v ’ v {t) := (C/ U {t),..., Cl/ U (t)), where C/ U {t) := ^ Y] for alii £ {1,..., m}, 

3 = 1 

where {{Yf : n > 1} : i £ {1,... ,m}} are m independent sequences of i.i.d. positive integer-valued 
random variables, and independent of {N' 1)V (t) : t > 0} as in Definition \l.ll 

Obviously the process {C v,u (t) : t > 0} in Definition 11.21 coincides with {N v,u (t) : t > 0} in 
Definition 11.11 when all the random variables {{1^ : n > 1} : i £ {1,... ,m}} are equal to 1; see 
also Remark [LT] below. In view of what follows it is useful to introduce the following notation. We 
start with the state probabilities {{q/ v {t) : k > 0} : t > 0} defined by 

q/ v (t) := P(C v,u (t) = k ) for all integer k\,...,k m >0, (3) 
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the probability mass functions 


qtj := P(Yf = j ) for all integer j > 1 (i € {1,..., m} and n > 1) 


and the probability generating functions 


Gi(u) := ^2 uJ Vj (* G i 1 ’ ■ ■ •> m l) 


i>o 


and 

We remark that 

G^iu-t) := E 




fci { 


■u- q r(t). 


k> 0 


U 


CiiAWm C m (A r '(C , '(t))) 

Urn 


= E 


E 


u 


ClM . . . u Cm(r) 


and E 


Ci(r) Cm(r) 
u x ■■■Um 


J r=A r l(L v {t)) 

= e^*=i^*( G i( u i)- 1 ) T1 ; thus, by taking into account (3.8) in [2J, we get 

(u;f) = E„,i ^ Ai(l - G*(^))j . 

As a particular case we can consider the probability generating functions 

C'-fet) :=X>l I ---“5n , rw 


( 4 ) 


and we have 


G^iuit) =E 


fc >0 


3 Er=i 


= ^,i (- (Ai(i-■«*)) e 


( 5 ) 


s, Z— 1 


note that both (JU) and ([5]) can be seen as a generalization of (3.20) in [2]. Finally we consider the 
probability mass functions concerning convolutions, i.e. 


(<f )* := P(Y\ + • • • + Y% = j) for all j > 1 (i G {1,..., m} and n > 1). 

We remark that, since the random variables {{Y„ : n > 1} : i G {1,... , m}} are positive, we have 

(cf )*° = l {j=0} ; if j < h, then ( q l )* h = 0. 

Remark 1.1 Obviously the state probabilities {{q]f v (t) : k > 0} : t > 0} reduce to {{pl’ v (i) : k > 
0 } : t > 0} when we have q ) := l{ J= i} for all i G {1 ,..., m}. 

A further generalization of the process {N v,u (t) : t > 0} in Definition 11.11 is the multivariate 
version of the process in p!6| . In view of this we recall that, given a nondecreasing Levy process 

(subordinator) {1H7 (t) : t > 0} associated with the Bernstein function /, we have 


E 




— e tf(n) ^£ or a £j n,t> 0); 


moreover we have the following integral representation 

/‘OO 

/(aO = (1 - e~ tir )p f {dr) (for all p > 0 ), 

J o 

where pf is the Levy measure associated with / (we also recall that pp is a nonnegative measure 
concentrated on (0, oo) such that J 0 °°(r A l)pf(dr) < oo). 
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Definition 1.3 Let us consider the processes in Definition [ITT] and an independent subordinator 
{!Kf(t) : t > 0} associated with a Bernstein function f. Then let { NL v (t) : t > 0} be defined by 




Remark 1.2 If {TCf (t) : t > 0} is the stable subordinator {APifb) : t > 0} cited above, we have (see 
e.g. Example 1.3.18 in m 


f(p) := fjfi, or equivalently Pf(dr) 


r(i - v) 


i 

r ri+l 


1(0,oo ){r)dr. 


Obviously in this case {N*' v (ft ) : t > 0} in Definition \1.3\ coincides with { N v,lJ (t ) : t > 0} in 
Definition \l . 1[ 


In what follows all the items concerning the process {Nf’ u (t) : t > 0} will be a modification of 
the ones for {IV 77 ’ 1 '(f) : t > 0} in Definition 11.11 with / in place of ry. thus, for instance, we set 

p{. ,l/ (t) := P(N^ u (t) = k) for all integer k \,..., k m > 0 ( 6 ) 

and 

Gt’-iu (7) 

fc >0 

We conclude with the outline of the paper. We start with some preliminaries in Section [2l The 
results are presented in Section OJ which is divided in two parts: 

1. the results for the processes in Definitions 11.11 and 11.21 

2. the results for the process in Definition 11,31 

Some examples of fractional compound Poisson processes and the generalization of a result in [3J 
for the fractional Polya-Aeppli process are presented in Section |4l 


2 Preliminaries 


We start with some useful special functions. We start with the generalized Mittag-Leffler function 
which is defined by 


E l/ X ) ■= Y. 

j> o 


j!P (aj + py 


(see e.g. (1.9.1) in [ 8 j) where 


( 7 )0' ) 


7(7 + 1) ••• (7 + j - 1) if j > 1 

1 if j = 0, 


is the rising factorial, also called Pochhammer symbol (see e.g. (1.5.5) in [ 8 ]). Note that we have 
E x a p, i.e. a with 7 = 1 , coincides with E a ^ in (fT|). 

We also recall the Fox-Wright function (see e.g. (1.11.14) in ( 8 ]) defined by 


fly ( a l) Oifi) ■ ■ ■ (flp, CXp) 

p q [ (6 1 ,/3 1 )...(6 (? ,/3 9 ) 


( \ . = se nLi + a hj) z j 

huUi m+Mfl’ 


under the convergence condition 


q p 

^ ^ fik ^ ^ &h E 1 

k= 1 h= 1 


( 8 ) 

(9) 
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(see e.g. (1.11.15) in 0). 

We conclude this section with the definitions of two fractional derivatives and of a fractional 
difference operator. Firstly we consider the (left-sided) Caputo fractional derivative of order v £ 
(0,1], i.e. c D” + in (2.4.17) in [ 8 ] with a = 0: 


C D v 0+ f{t) := 


r(i^y/o JTsrTsf( s ) ds if ^€( 0 , 1 ) 

if v = 1 - 


( 10 ) 


We also consider the (left-sided) Riemann-Liouville fractional derivative of order u > 1 (see 

e.g. (2.2.4) in 0) defined by 


d v 

d(-ty 


m ■■= 


(- m) ' n It°° if v is not inte § er and m - 1 < v < m 

(-1 if v is integer. 


Moreover, for r/ £ (0,1], we consider the (fractional) difference operator (I — B) v in [15]. More 
precisely / is the identity operator, B is the backward shift operator defined by 


Bf(k) = f(k-1) 


( 12 ) 


and, if we consider the Newton’s generalized binomial theorem for operators, we have 


(i- B y = Y^(-iy 

3> 0 



IP. 


3 Results 

In general we show that the state probabilities (and the probability generating functions) solve 
suitable fractional differential equations and we provide some explicit expressions. In order to have 
a simpler presentation of the results, throughout this paper we always set 

m 

«(A) 

i=l 

(also in the next Section HD, where A = (Ai,. .., A m ). Moreover let {Bi : i £ (1,..., m}} be the 
operators defined by 

Bif{ki ,..., k m ) = f(k i, — 1,..., km); (13) 

these operators play the role of the operator B in (11211 for the case m = 1. 

3.1 Results for the processes in Definitions 11.11 and 11.21 

The first result shows that the state probabilities {{p'l' v {t) : k > 0} : t > 0} in ((2D solve fractional 
differential equations, and we consider the fractional derivative in ( 11011 . 

Proposition 3.1 For rj, v £ (0,1], the state probabilities {{p^ v {t) : k > 0} : t > 0} in ((2} solve 
the following fractional differential equation: 

f {' - £ i^ £i )%r(*) 

1 pT(t) = i( 4 -o>. 
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Proof. Firstly, by ©, we have 


C DZ + G^(u; t) = - (££1 Aj(l - m))" G v,u (u] t ) 
G v ' u (u; 0) = 1 


by (2.4.58) in f8j, and therefore 
[ °D? 

\ G r >’ v (u;0) = 1 . 


o + G^(u;t) = -(«(A)r (l - 


(14) 


From now on we concentrate the attention on the first equation only (the second one concerning 
the case t = 0 trivially holds). Then, if we use the symbol £ r m es f° r H ie sum over 
ri,..., r m > 0 such that r\ + • • • + r m = j, we have 


1 - 


«(A) 


j> o 


E (-F 


/ sr^m \ 

I £?:=i Ai'U, 


.7/ V »(A) 


y {*) ( ~ 1)J y 


i>o 


J 




rAf " 1 •••Air -u? •••«; 


j7(*(A)V ri 1 m 1 


Thus 


c os + ff'--teo = -(*(A))"E())(^ E 

i>0 W V ri,...,r m eS. 


n! •• - r m ! 


A?-AJ 


Y J 4'+ r '-’^-‘pr(t) 


k> 0 


where, for the last factor in the right hand side, we have 

£ »»■+'■ ■ ■ ■ = E ”1‘ " • “m p^ r w- 


fc >0 


k>r 


Then (in the next equality we should have r\ < k\,..., r m < k m , but this restriction can be 
neglected) 


m; i) = - («(A))" ■ ■ ■ <*• j] (A 

fc>0 j >0 


/? A (-*7 


E 

ri,...,r m e§j 


n! • ■ ■ r m ! 


j7 7(A)7 

A?---A^r r w- 


We conclude the proof noting that, since 


E 

ri,...,r m eSj 


r\\ 


A? • • • KnPk-r^) = ( J2 XiBi ) pY®' 


K 1=1 


where B i,..., B m are the shift operators in (1 13 [i . we have 


C 


D^G^(u-t) = -(s(\)ryu k 1 


k l . . , n.kn 


k> 0 


Eli 

3> 0 


v\ (-17 


jJ (s(A)7 


E A ‘ B d prw 




(»(A))"E“f 


fcl ( / - 


k>0 


£ Zi 

s(A) 
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which yields the desired equation. □ 


The second result concerns the state probabilities of the fractional compound Poisson process, 
i.e. : k > 0} : t > 0 } in (J5J). More precisely we mean {{ql' u (t) : k > 0 } : t > 0 } (time 

fractional case) and {{ql’ u (t) : k > 0} : t > 0} (space fractional case). We show that they solve two 
fractional differential equations: the first one is a generalization of Proposition 13.11 with q = 1; in 
the second one we have the fractional derivative (HD. 


Proposition 3.2 For v € (0,1], the state probabilities {{q^it) : k > 0} : t > 0} in ([3]) solve the 
following fractional differential equations: 


C D v 0+ q]f u (t) = -s{\)q]f v 
Qk U ( 0 ) = l { fc=o } - 


(t)+Er=i 




For rj £ (0,1], the state probabilities {{(/^(t) : A: > 0} : f > 0} in (T3|) solve the following fractional 
differential equations: 


3^4' 1 M = s(AW' 1 (*)-£E,i 
qf(0) = 


\ , n l n ,A 

A * 2->ji = 1 ki-ji, 


V’ 1 




it) 


Proof. Firstly, by (ID), we have 

f c D- + G 1 6 l '{u-t) = -ET =1 Ai(l -GiiuiVG^fat) 

1 G^(u- 0) = 1 

by (2.4.58) in [ 8 j and 

( = Eii Ml - 0,M)Op( u: t) 

\ cj'fe 0 ) = 1 

by (2.2.15) in [ 8 |. In both cases the second equation (concerning the case t = 0) is trivial, and 
therefore we concentrate the attention on the first equation. So, if we compare the equations above 
and the ones in the statement of the proposition, we have to check that 


m ki 


E Ai(l - Gi(ui))G^(u; t) = J2 4 1 ■■■«£?( sQdqi v (t) + E Ai E ^kZ-M-F-.kJ^ 

k> 0 \ i= 1 ji= 1 


i =1 


and 


X>(1 -G i (u,))Gf{ m t) = Y J u\ 


k\ gi^ rn 
a m 


(A )Qk\t) - E E %i*ku..,ki-3i,.,kn 
i= 1 ji =1 


it) ; 


i= 1 fc> 0 

moreover, after some easy manipulations, the above equalities are equivalent to 


m m ki 

E AjGi(u i )G ( 4 (u; t) = E u ~i ' • • u m n E A * E q)i < tk 1 ,..;ki-ji,-,K 

k> 0 i= 1 ji= 1 


it) 


%— 1 


and 


E AiG , i(u i )Gj 1 (u;t) =E 


,&1 


1=1 


k> 0 


E A < 

2—1 


E 

Ji=l 


~i r?,l 


it), 
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respectively. In the first case we have 


Y^\Gi(ui)G 1 ^{yr,t) Ai u ?4 5Z “l 1 " ' v fc<&''( t ) 

i= 1 3 i>l fc >0 


t= 1 


i= 1 ji>l k>0 


and the desired equality holds because the sums and the factors in the last expression can be re¬ 
arranged in a different order and . fcm W = 0 when j, > hi- The other case can treated 

in the same way (we have to consider G ’^ 1 and {(g)?’ 1 ^) : k > 0 } : t > 0 } in place of G]f u and 


As a special case we give a version of the equations in Proposition 13.21 for the state probabilities 
{ {p r j/ / (t) : k > 0} : t > 0} in (J2J) for the multivariate fractional Poisson process in Definition 
11.11 The first equation meets Proposition 13.11 with rj = 1; the second equation with 7 / = 1 meets 
Proposition 13.II with 77 = v = 1 (i.e. for the non-fractional case). 


Corollary 3.3 For v G (0,1], the state probabilities {{p]f v (t) : k > 0} : t > 0} in ([2]) solve the 
following fractional differential equations: 


c osy/(t) =-«( a + e ” , yi 

= 1 {k=0)- 


,ki 1 5 • • • 5 kfYi 


(t) 


For rj G (0,1], the state probabilities {{p^it) : k > 0} : t > 0} in (|2|) solve the following fractional 
differential equations: 


Pfcl 0 ) = !{fc=o}- 


Proof. It is an immediate consequence of Proposition 13.21 and Remark ll.il □ 


Now we give some expressions of the state probabilities {{Pp.' 1 '(t) : k > 0} : t > 0} in (]21) . We 
start with an implicit expression which generalizes (3.19) in [2] ("note that we use the notation d\ i 
in place of ^-). The most explicit formulas are given in Proposition 13.51 


Proposition 3.4 Let 7/, v G (0,1] be arbitrarily fixed. Then, for all integer k \,..., k m > 0, we have 

m 

Proof. By construction we have 


1=1 


P Y{t) = E 


E 

n{ (A :fW 



.41 kf. J 

Li= 1 v x 

z=A r > (£ v (t)) 


1 


fei! • • • k m \ 


-E 


m 

no-w *} 1 


-s(X)A-n (£"(t)) 


_i =1 


then we can conclude by following the same lines of (3.19) in [2], where we take into account that 
E [ e -RA)^(£-(t))j = EvX (-( S (A))^) by (3.8) in [2]. □ 


Proposition 3.5 Let 77 , v G (0,1] be arbitrarily fixed. Then, for all integer k \,..., k rn > 0, we have 

(_i)fci+-+fc m ^ (-(s(A )m T(gr + 1) 

- (a( A)) fe i+-+ fe - kfi- ■ ■ k m \ ^ T(vr + 1) r(r)r ~ (h + ■ ■ ■ + k m ) + if 1 j 
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or equivalently 

X'l 1 • • • A 


PW = 


(- 1 ) 


fclH- 


(s(A)) fc i+-+ fc - k\\ - ■ ■ k m ! 


2^2 


(M) 


( 1 , 1 ) 


(1,^) (1 — {ki H- \-k m ),v) 


i-(sm v n 

(16) 


Proof. The equality (fT6l) follows from (fT5l) . In fact, by taking into account ([8]), it suffices to multiply 
the terms of the series in the right hand side of (11511 by = 1 (note that the convergence 

condition ([9|) holds because u + rj — (ji + 1) > —1). So from now on we can concentrate the attention 
on the equality (fT5l) only. 

Firstly we have 


pT(t) =P (= i}n|x>rw = 

/ m m \ / m m \ 

=p N v,u (t) = k\j_i^(t) = h • p E = 


(17) 


2 = 1 


2—1 


v 2— 1 


2=1 


We start with the conditional probability in (117[1 : then we have 

P (N r,,L '(t) = k) 


p[N^{t) = k\Y J Nr(.t) = Y J ki 


2=1 


2=1 




and, if we consider the conditional distributions given A ri (L u (t)), we get 


P [N^(t) = k\Y^N^(t) = J2h 


E 


nr=i 


(\ir) k i —Xif 
ki! e 


r=A’i(£ , '(t)) 


2=1 


2=1 


E 


(s(A)r)gfal (A)r 


(ESTfcOl 


r=A r >(£ ,/ (t)) 


{k\ + • • • + k m )\ A^ 1 • • • A^ 


k\\---k m \ (s(A)) fcl+ "' +fc " t 

after some computations (there is a factor equal to 1 given by E (A'' 1 (L u k i e - s (P) AV {*)) 

divided by itself). For the second factor in (11711 we consider again the conditional distributions 
given A r} {L v (t)) and we have 


p =e 


< 2=1 


2=1 




< 2=1 


2=1 


r=A'n(L u (t)) 


then we get 

/ m m 

p 


(-i) 


=E 


&1 H-h k-n 


(s( A)r)E^ 7 

(E)=i h)\ ' 


r=A , n(L u (t)) 


< 2=1 


2=1 


(k\ + • • • + k m )\ 


E 

r>0 


(-(«(a )rn 


T(j]r + 1) 


T{vr + 1) r(r/r — {k\ H-h k m ) + 1) 


by taking into account the known formula for the case m = 1 (see (3.24) in [2] where the formula 
is given in terms a binomial coefficient and there is a typo; see also (1.8) in [15]). Finally (11511 can 
be easily checked. □ 
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Here we present some remarks on Proposition 13.51 Firstly (1151) with m = 1 meets known 
formulas in the literature (see e.g. (1.8) in [15]^. Moreover, for v = 1, we have 


\k\ \k. 

pf(t)= A1 '" A " 


(- 1 ) 


kl~\ - \~kn 


(s( X)) kl+ '" +km ki\■ ■ ■ k m l 


E 

r>0 


(-(^(A))^) r 


r (rjr + 1) 


r\ 


r(r/r — (k\ H-b k m ) + 1) 


and 


Pf(t) = 


\ ki \k 
A, • • • A r 


(- 1 ) 


kl~\ - \-kn 


(s(A)) fcl+ '"+ fcm k\\ - ■ ■ k m ! 


i^i 


(M) 

(1 - (h H-1- k m ),T]) 




both formulas reduce to the ones in Theorem 2.2 in m concerning the case m = 1. Finally, for 
rj = 1, (fT5l) reads 


\ki . . 

\km ( _1 'ifclH-hfc 

_ A 1 A m V V 

Pk \ Z ) ~ 


E 


(-»(A)C)’ 


T\ 


(s(X)) kl+ '"+ km h\---k m \ r(i/r + l) (r — (fci + • • • + A™))! 

P K, 1 H-hfcrr?. 

(because the summands with r < fci + • • • + k m are equal to zero), and therefore 

(_i)fci+-+fem (-s(A)^) r+fc i+'" +fem {r + k\ -\ - \-k m )\ 


P?(t) = 


E 


(s(A)) fcl +-+ fcm Aq! • • • fe m ! ^ T(z/r + z/(Aq 4-bfc m ) + l) 


Y i 


(fcl + • • • + /c m )! 


\All \km j.i 

Ai • • • A™ • l 


'(fcl+-"+fcm) . (^1 H + km + l) (l) • { — s(X)t U ) r 

^ d-r^r + ^l-f A: m ) + 1) 

^ (fcl + • • • + fc m ,)! ^ \fc m 4.v{k\-\ -h km) 171(^1 4 -bfcm) + l 7 n l\\4-V\. 

-—- A 1 '" A ™ ■^( fc i+...+fe m )+il- s (A)t ), 

the last expression meets (2.5) in [5] concerning the case m = 1. 

In the next Proposition 13.61 we compute the covariance 


Cov #(1),^ 


1 , 1 / 


:= E 


N^(t)N^(t) 


-E 


lY^(t) E IV^(t) (for j, he m}); 




note that we take ry = 1 otherwise the covariance would not be finite. In what follows we refer to 

1 


m--=- 1 


v \T(2n) vT 2 (v) 


(18) 


where, as shown in [3] (Subsection 3.1), Z(y) > 0 for v G (0,1] and Z(v) = 0 if and only if 
v = 1. The codifference t(Xi,X 2 ) is studied in the literature (see e.g. (1-7) in [9j) when the 
random variables X\ and X 2 have infinite variance and it is known that it reduces to Cov(Ad, X 2 ) 
when ( X\,X 2 ) forms a Gaussian vector (see the displayed equality just after (1.7) in [9]). So in 
Proposition 13.61 we also compute the codifference 


: = logE 




— logE 


e m rw 


— logE 


0 -iNZ >v (tj 


(for j,h G {l,...,m}), 


where i is the imaginary unit. 

Proposition 3.6 Let r),is G (0,1] be arbitrarily fixed. Then, for j, h G {1, ..., m}, we have: 


Co w[N]' v {t),N L h 


1,IS 


^{j=h} ’ 
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where Z{u) is as in (USD; 

r (iVf-(t),JVj'-(t)) =l Wh) ■ log -«*) + A»(l - c- , ))”n 

- logE^HWl - e‘))V) - log^,,(-(Ak(l -e-^V), 

where i is the imaginary unit. 


Proof. Firstly it is useful to recall the following formulas: 

X k t l 


E 




T(n+1) 


(for all k G {1,... , to}) 


(see e.g. (2.7) in [1]); 


E 


e iuN]l' v {t) 


= E U) i(—(Afc(l — e m )) ,? t^) (for all u € M and k € { 1,... , to}) 


( 19 ) 


( 20 ) 


which can be obtained by adapting the computations in m for the generating functions. 

We start with the case j = h. The formula for the covariance holds noting that Cov(lV J 1 ’ I/ (t), Nj’ u (t)) 


Var 


Nj’ v (t) and by taking into account (2.8) in [5]. The formula for the codifference holds noting 


that E 




= 1 and by taking into account (1201) . 


We conclude with the case j ^ h. Firstly we have 


E 




= E 


E[-^j 1,1 ( s )]®[lV^’ 1 (s)] 


s=L u (t) 


j 


= XjXh S f^(t){s)ds 


and, since 


f 


U\ . juk 

s k fc»(t)(s)ds = (for all k > 0) 


T(vk + 1) 


by combining (2.4) and (2.7) in [17| . we have 


E 




— ^h 


2 1 


2u 


then, by taking into account (fT9lh we obtain 
Covf Nl’ 1 '(t),Ni' v (tj) =XjX h 


2 1 


2v 


r( 2 ^ + i)’ 


Xjt h 


\ h t v 


' y ' n T(2u + l) T(u + l) T{u + 1) 

2 1 


—\j 


2 v 


T(2v + 1) r 2 (i/ + l), 

X]Xhf2v (2^(2^) “ }^r>)) = w 2 "^) 

and the formula for the covariance is proved. Furthermore, since we have 


E 


=E 

E 

IN}’ 1 (s) 
e 0 v ’ 

E 

e ~ iN h 1 ( s ') 








- 

s=A r i (A" (t)) 


=E 


0 \js{e l -l)+\ h s(e l -l) 


= E v>1 (-( X j (l-e i ) + X h (l-e- i )) 1 >1?), 


s=A r > 

the formula for the codifference can be easily obtained by taking into account (j20lh □ 
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It is known that {C v,1 (t) : t > 0} and {1V 7?,1 (£) : t > 0} are Levy processes and, moreover, when 
rj = 1 their Levy measures p l c and p]^ are defined by 


Pc (A l x ■■■ x A m ) — 'y ' XpffAf) 


2—1 


and 


Pn{A i x ••• x kf m ) — Ail{i ej4 . } . 


( 21 ) 


( 22 ) 


2—1 


In the next proposition we present the Levy measures p n c and p)y when p € (0,1). 

Proposition 3.7 Let p € (0,1) be arbitrarily fixed. Then the Levy measure p n c of {C 71,1 (t) :t> 0} 
is defined by 


P(j{.A i x ••• x ^4 m ) — 


E / n E « 




72 j >0 




, (A^) n “ 


m ! 


Moreover the Levy measure p v N of {N 71,1 (t) : t > 0} is defined by 


i x ••• X Al m ) — 


E 




r(/ci H-b k m - p) f Aj* 


3 -s(a)z 

^ +1 


-da. (23) 


r(1 “ ^ fcl ( s (^)) fcl+ "' +fem_?? i=1 


n i fe-i ■ 1 


J ' HLeA;} 


(24) 


Proof. Firstly, by (30.8) in )2T and the Levy measure pj for the stable subordinator {A u (t) : 2 > 0} 
in Remark 11.21 we have 


pl(A 1 x ■ ■ • x Am) = J2 / I] E 


ky 0 " ^ i= 1 I Ui>0 


li (Aj-2) 


'fc; 


; AjZ 


re,;! 


• 1 


{ki&Ai} 


r(i - r?) z^ +l 


dz. 


Then we easily get (1231) with some manipulations. Finally, as far as (1241) is concerned, we have to 
consider (1231) with g*- := l{j=i} for all i £ {1,... , m}; therefore we have (rf)£ n * = l{ k . =ni } and we 
obtain 


x ••• x A m ) — 


V 


r(l-ry) 


/»oo m 

e/ n 

ky 0 JU i=l 


(A*z 


. >y\ki 


ki\ ' 1{fcieAi} / Z V+1 


D —s(X)z 


-dz 


” l x k 






2—1 


which yields (1241) . □ 


We remark that in (1231) meets (1211) . In fact, if we set = 1, we have a non-null 

contribution if and only if (n ±,..., re m ) belongs to the set {(1,0,... , 0),..., (0,... , 0,1)}; thus 
yields 

1 POO m 

Pc(A 1 x X ^4 m ) = 

A i =1 fci>l 


1 r(i-i) 


r(i - 1 ) (s(A))° 


E A < E Rw*>} = E a <«a>. 


2—1 ki> 1 


2—1 


Similarly pjy in (1241) meets ()22f) . In fact we have a non-null contribution if and only if (k ±,... , k m ) 
belongs to the set {(1, 0,..., 0),..., (0,..., 0,1)}, and (l2ll) yields 


A(AiX...x^) = F j i L_.fffigj».W 


(leAR — 

2—1 
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3.2 Results for the process in Definition 11.31 

Here we give a multivariate version of Theorem 2.1 and Remarks 2.3 and Remark 2.5 in [16j . In 
particular we recover those results and remarks by setting m = 1. In view of what follows we 
consider the analogue of (1.1) in [16], i.e. 


P(N f ’\t + dt) -N f ’\t) = k) 


/ 0 °°(n,=i (A k < 1 e Xir )pf(dr)dt + o(dt) for k >- 0 

1 - /o°°( 112=1 e- Xir )p f (dr)dt + o(dt ) for k = 0 

ni=i XT ' fo° r ^= 1 ki e~ s ^ r pf(dr)dt + o(dt ) for k >- 0 
1 — / 0 °° e _s (-) r pf(dr)dt + o(dt ) for k = 0 


and we consider the function f m defined by 


/™<A; u):= 

j>0i=l 




ji ! 


)Pf(dr ); 


in particular we have 

POO 

fm( A;0)= / (l-e- s ^ r )p f (dv) = f(s(X)) 

■Jo 

for u = 0, and 

POO 

fi(Xr,u 1 )= / (1 - e~ Xir+Xluir )p f (dr) = /(Ai(l - «i)) 

Jo 

for the univariate case m = 1. 


Proposition 3.8 Let f be a Bernstein function. Then we have the following results. 

(i) The state probabilities {{p{f l (t) : k > 0} : t > 0} in ((6]) solve the following fractional differential 
equation: 


TtPk 1 ^) = Eg^jKkPk-jiV W?=i jfr r r^ ji e rs ^p f {dr) - /(s(A))pf 1 (t) 


Pk 1 ( t ) = 

(ii) The probability generating functions {G-JJf^t) : t > 0} in © solve the following fractional 
differential equation 


tit 


G^(u- 0 ) = 1 , 

and therefore we have 


Gf' 1 (u;t) = -f m (X] ujG^iu;1) 


Proof. We start with the proof of (i). The initial condition trivially holds. Then, since {N-fJ^t) : 
t > 0} has independent increments, by taking into account the distribution of the jumps given 
above we have 


p{’ 1 (t + dt)= ^2 + dt) - = k — j) 

o <l<k 

= J2 Pf'^lj TI j.y e ~ Xir )pf( dr ) dt + o{dt)\ 
0 <i<k V 1 ' 0 i=1 : 1 dtp / 

+ p{’ 1 (t) ^1 - J e -s( - )r p f (dr)dt + o(dt)^j , 
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and therefore (we consider a suitable change of summation indices in the last equality) 

/ m \ ki —ji /‘oo 

pi’ 1 (t + dt) - pi’ 1 (t) = ^2 Pj :1 {t) ( II (U l _ ■ \| / r^i^ l( - ki ~ ji ^e~ s( ^ )r pf(dr)dt + o(dt) 

" \<=i 1 * • 7ij - 1/0 


0</ -< fc 


- Pfc (*) (f(s{X))dt + o(di)) 

/ m \ji /‘OO X 

y: Pfc-j W I n “u" / r'£' i = lji e- s( 'X) r Pf(dr)dt + o(dt.) 


0<j <k 

JA 


\i 1 


Jd Jo 


-Pk (t)(f(s(X))dt + o(dt)). 

We conclude dividing by dt and taking the limit as dt goes to zero. 

Now the proof of (ii). The initial condition trivially holds. Then, if we take into account the 
differential equation obtained for the proof of (i), after some manipulations we get 


j t Gf 4 («; t)=J2 “ 1 1 ■■■ uk m J t Pk 4 (*) 


/c>0 


( m \ Ji /*00 

2 / r ^= l3ie ~ rs{ ^pf( dr ) - f{s{x))p f /{t) 

EIK ( E *Wnf rr^n^.W pf{d r)\ ; 

k>0 i=l \0-<j<fc i=l </0 / 


= -/(^(A))G / - 1 (u;i) + 


moreover, if we rearrange the summands in a different order, we obtain 


d_ 

dt 


m ( m \ 3i r°° 

G fd (u]t) = - /{siX^G^^t) + J2J2H u i i U&WII-T / ^ EI = lJi e- rs(A) p/(^) 

r*oo 




/*00 ( \ \ 7' 1 ^ 

= -/(s(i))0 , ’ 1 (a;*) + E / *-”“[1^o/WEII 

/M)-' 0 i=l ■'*' fc>2 i=l 

/MA)) + £ r e-"® n G/ ’ lfe t); 

2>-o^° j=i ■'*' ) 

finally we can check that (in the first equality we take into account the integral representation of 

/> 

r e --a)jj(Wh w(dr) ] G /,!(„.() 

W° J>0^° i=l / 


_d 

dt 


3 fYl / \ \ d • 

(1 - e-<a ■ V n ^^)p f (dr) ] G"(„; t) 


j>0 i= 1 


Jd 


= - fm (A; it) G f ’ 1 (u; t), 


and this completes the proof. □ 

Remark 3.1 The equation in Proposition \3.8f i) can alternatively he written as 







where B_ = {B \,..., B m ). In fact we have 


= - /”(i - e-®- 

7 ° j>0 i=l 


°° —s(X)r \ f,lu\ 

e z_.ll— 71 —p/( 7r K w 


j>-0 i=l 


-/(s(A))Pfc 1 (i)+ [ 

Jo 

m \ji POO 

II “V / r^^ Ji e~ s ^ )r p f (dr) - /(•s(A))p{’ 1 (i) 

>0 " i=i Ji ' • y ° 

m \ji POO 

J2 Pfcwn-T / r^= lJi e~ rs< '-^p f (dr) — /(s(A))p{’ 1 (t) 
-<7<Jfc " i=l Ji ' - 70 


0-<J < fc 

Remark 3.2 If we follow the same lines of Remark 2.5 in /or z/ € (0,1) i/ie state probabilities 
{{Pk’ U (t) : A A Q} : t > 0} m ([6]) sofee i/ie fractional differential equation 

c 0^4"it) = - Hs(\))p r P(t) 

pfrit) = 1 ( 4 - 0 ). 


or equivalently 


° DZ + pfc V (t) = - fm{kB)p f k'(t) 
pl’ U (t) = !{fc=0}- 


(25) 


Moreover the probability generating functions {G pu f;t) : t > 0} in ED solve the fractional differ¬ 
ential equation 

C D" 0+ Gf’"(u-1 ) = -/ m (A; u)G^(u; t ) 

g^(m;0) = i, 


(26) 


and therefore we have Gf ,u (u-,t) = E Ut i(—t u f m (X; 1 u)). 

In particular, if we consider the Bernstein function f for the stable subordinator {A v (t) : t > 0} 
and the corresponding Levy measure pf (see Remark \1.2\) . we have 


fr 


roo 

»(A;«) = / (1 

J 0 


_ p -s(A)r 


eii 


(A iUir) ji 77 


i>o 


f = , J'J r(l-r/) r 


7/+1 


dr 


=(«(A))" - 


=( s (A)r + 


v 


-7 ? r(-7?) 


oVn 7—1 ^0 


r-EIii Ji —»7—1 e —s(A)r- 


jVO i=l 


=(»(A))’ ^1 + 

-ma^E 




r (-r?) 

r (EiAl A - ? ?) TT 1 


n/ 


i>0 jLl^WA). 

moreover, if we use the symbol j m e§ h f or th e sum over all j±,... ,j m > 0 such that j\ -\ -+ 
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jm = h (as in the proof of Proposition \3.1\) . we obtain 


frn&u)=(sm V Yl 


r {h-rj) 


e m 


hi f XiUi\ 7i 


T(-i])h\ ^ 1L ji'-\s( A) 


=(s 


(e a ™ 


h>0 


r( ??)/i! \^a(A) 


/ m . \ ^ 

= MA))” i-E 


Et 5 0>, 


([/or f/ie last equality see e.g. (15) in \23f with a = —77 — 1 and /3 = 0; in fact t and ( in that 
reference satisfy ( = t(l + £), and therefore ( = and 1 + £ = ; obviously here we consider 
ui,...,u m G [0,1] and therefore t = G [0,1],). Thus (1251) meets the equation in the 

statement of Proposition \3.1\ (with p r () u {t) in place of pjj 1 ' (t)) and, similarly, (1261) meets (TH1) (with 
G v,,y (u]t) in place of G^ ,u (u-,t)). 


4 Examples of fractional compound Poisson processes 


In this section we study the multivariate fractional version of well-known counting processes which 
can be obtained as a particular multivariate space-time fractional compound Poisson process 
{C n,u (t) : t > 0} as in Definition 11.21 In particular the univariate processes (i.e. the case m = 1) 
has been studied in |3j (Section 4). For each example we specify the probability mass functions 
m-- j > 1} : i € {1,..., m}} and the values Ai,..., A m ; we remark that the values Ai,... , \ m in 
Example 14.11 can be chosen without any restriction. 

Example 4.1 (Multivariate fractional Polya-Aeppli process) We set 

q) := (1 - 

for some di,... ,a m G (0,1]; in particular, if dti = 1, we have = NT v (t). We recall that in 

some references the case m = 1 is presented with p in place of 1 — a; see e.g. (1.3) in Ill- 

Example 4.2 (Multivariate fractional Poisson inverse Gaussian process) We set 

(3-3/2) ( 2ft V 

fj := - J P/W ani := f ((1 + ~ 0 

( 1 ) 7 _ 1 Pi V J 

\2/3i+l) 

for some (3i,fi x ,..., /3 m , fi m > 0. 


Example 4.3 (Multivariate fractional Negative Binomial process) We set 

(1 - di ) 7 


q) ■= 


and \i := — log d* 


3 ' j log dj 

for some di,.... d m G (0,1). 

We also present an extension of Proposition 2 in [3] concerning Example 14.11 
Proposition 4.1 Assume to have the situation in Example \). /[ Then: for v G (0,1], 


c o 5 + «P(*) - £,”1(1 - «.) • .JD 


= -■(AtepM + ES.(AA + »(A)(1 - a,)K:...,h-i .*.<*> 


1m 


- E”i(i - a.)EZU,**ES„(i - .*„(*) 


l.i/ 


9^(0) = i { fc=o } ; 
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for i] G (0,1], 


^(*) - 5Xi(i - *0 • srU!W) 

= *(A)?f 1 (t) - £r=i(A^ + »(A)(1 - «*))<..., 

+ E’;=i(l - a<) Y2=\ y h+i x h E^=i(l - W 

gf(0) = l { fc=o } - 

Proof. The initial conditions trivially holds. We start with the proof of the first equation in the 
statement. By the first equation in Proposition 13.21 we have 

m 

C n*' 




i= 1 


m kf,, 

sQdqi v (t )+E A/ > E t 1 _ 

/i=i jh=i 


Et 1 _<5 *) 


2=1 


kh jhi---)kr 


kh 


,(*) 


h=l,h^i jh=l 


+\i 5^(i - di) Ji 1 Oig^... ifci _i_ iJ - 4) ... )fcm (t) 

h=i 

Moreover, if we split in two parts the sum ^^ =1 (1 — k h -j h fc m (A) i n the right hand 

side, i.e. the summand with j/i = 1 and the other summands with € {2,..., after some 
computations we get 


C A)+<7fc’ (t) E(-*- “j) • 


2=1 


m m kh 

s (A)i’"w+E + E A)i D 1 - 

h=l jh=2 


/i=l 


+ 5E«(A)(1 - 




E^ 1 ~ “A E A h 5E - 


2=1 

m m 

- &i) 

i=l h=l,h^i jh =1 

m ki 

E Ai E(f - .fc™ (*) ■ 

2=1 J*=l 


kfi jhi'"ik‘r 




Finally, after some other computations (in particular we put together two sums and we consider 
ji € {2,... , ki + 1} in place of ;j t € {!,..., A:*} in the last sum, with a suitable modification of the 
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summands), we have 


c m + ii' , it) - ][> - a,) • c Dz +q i^ ti _ 1 . k jt) 


i =1 

m 

- s(\)ql u (t) + i&i + a(A)(l - w))q 1 k*..., ki -i,..., km (t) 

i= 1 

m k h 

+ 12 Xh Jl^- & *Y h ~ 1& h *k?,.,k h -3 h ,.,km® 

h =1 jfc=2 

m m kh 

- 5^(1 - an) ^ A ^ ( X - ^) jh ~ la ^kZ., kh - jh ,.., kjt) 

i= 1 h=l,h^i jh = l 

m ki~\~ 1 

*=i h=2 


Then the first desired equation is checked because ^ fc k (t) = 0 and two sums can 

be canceled. The second desired equation can be obtained similarly; we have to consider the second 
equation in Proposition 13.21 (instead of the first one) and we have the same kind of computations 
with suitable changes of sign. □ 


Acknowledgements. We thank Bruno Toaldo and Federico Polito for some useful discussions. 
In particular Bruno Toaldo gave us several comments on the content of the reference m • The idea 
of studying the processes in this paper was inspired by the communication of Daniela Selch at the 
EAJ Conference in Vienna (September 10-12, 2014). 
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